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abstract 

The Maddenjulian oscillation (MJO) is the dominant mode of intraseasonal variability in the 
tropics. Despite its primary importance, a generally accepted theory that accounts for fundamental 
features of the MJO, including its propagation speed, planetary horizontal scale, multi-scale 
features, and quadrupolc structures, remains elusive. In this study, we use a shallow water model 
to simulate the MJO. In our model, convection is parameterized as a short-duration localized 
mass source, and is triggered when the layer thickness falls below a critical value. Radiation is 
parameterized as a steady uniform mass sink. Slowly eastward propagating (MJO-like) signals 
and red noise spectra are observed in our simulations. In the time-longitude domain, MJO-like 
signals with multi-scale structures are observed. In the Fourier domain, spectral peaks associated 
with the MJO-like signals are observed. In the longitude-latitude map view, quadrupole vortex 
structures associated with the MJO-like signals are observed. We propose that the simulated MJO 
signal is an interference pattern of westward and eastward inertia-gravity(WIG and EIG) waves. 
Its propagation speed is one half of the speed difference between the WIG and EIG waves. The 
horizontal scale of its large-scale envelope is determined by the bandwidth of the excited waves, and 
the bandwidth is controlled by number density of convection events. Our results suggest that the 
MJO perhaps is not a large-scale low-frequency wave, in which convection acts as quasi-equilibrium 
adjustment. Small-scale high-frequency waves might be important. 



1. Introduction 

The Madden- Julian Oscillation (MJO) is the dominant 
intraseasonal variability in the tropical atmosphere. It is 
a planetary scale, slowly eastward propagating (about 5 
m s~^) perturbation of both dynamical and thermodynam- 
ical fields. During an MJO event, a positive convection and 
rainfall anomaly develops in the western Indian Ocean and 
propagates to the western Pacific Ocean. Once the per- 
turbation reaches the date line, the perturbation is largely 
confined to dynamical fields. The associated planetary- 
scale wind structure is baroclinic, and is characterized by 



low level convergence and upper level divergence (Madden 



and Julian (19721; Madden and Juhan (1994); Hendon and 



Salby (1994)). More details can be found in comprehen- 



sive reviews, such as, Madden and Julian ( 2005 ) and Zhang 



(2005). In addition, within the MJO envelope, there are 



both westward and eastward moving fine scale structures 



(Nakazawa (19881, Hendon and Liebmann (1994)). 



The significance of the MJO for phenomena such as 
monsoon onset, ENSO, and rainfall patterns in the Tropics 
has been well established (Zhang (2005)). However, a suc- 



cessful MJO theory is elusive. For a historical and detailed 



review of theories, please refer to Zhang (2005) and Wang 



(2005). Here we will review three schools of theories. One 



school considers the MJO as a largescale unstable mode in 



the tropics, and it is often referred to as the moisture mode. 
The moisture mode arises from positive feedbacks between 
precipitation and the source of moist static energy, e.g., 
Neelin and Yu' (Il994|), iSobel et all (|2001|), iFuchs and Ray- 



mond] (|2002, 2005 ), 



Raymond (2007), 



Bretherton et al. (2005), Fuchs and 



Maloney (2009), Raymond and Fuchs 



( 2009[ ), [Andersen and Kuangj (|2012j). Another school still 
considers the MJO is a large-scale mode, but the major 
instability to maintain the planetary scale envelope hap- 
pens in the synoptic scales. Both observational and mod- 
eling studies suggest that high-frequency small-scale waves 



are important to the MJO, e.g.. 


Hendon and Liebmann 


(1994), Tung and Yanai (2002), 


Moncrieff (2004), 


Biello 


and Majda (2005 


), Khouider et al. (2012). Thus 


Majda 


and Stechmann ( 


2009 2011) emphasize the importance of 



small-scale waves within the MJO envelope, but they pa- 
rameterize the effect of the waves. In their model, it is only 
the amplitude of the wave activity envelope that is needed, 
not any of the details of the particular synoptic scale waves 
that make up the envelope. The third school considers 
the MJO as a wave packet of a certain type of equatorial 
waves, and the MJO propagates with the group velocity 
of the equatorial waves. Yang and Ingersoll ( 2011| ) hy- 
pothesize that the MJO is a mixed Rossby-gravity (MRG) 
wave packet that propagates with the MRG group velocity, 
and they test this hypothesis with both an idealized GCM 
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and with the Outgoing Longwave Radiation dataset. They 
force with a westward-moving heat source that lasts for 10 



or 20 days. At the same time, Solodoch et al. (2011) sug- 



gest that the MJO could be a forced response to the MRG 
wave group in the quasi-equilibrium (QE) convection and 
wind-induced surface heat exchange context. Neither the- 
ory considers inertia- gravity (IG) waves, which have fre- 
quencies greater than 0.5 cycles day~^. In the QE context, 
these high-frequency small-scale waves will be damped the 
fastest, and it is possible that the QE approach does not 
apply there. Therefore the role of high-frequency small- 
scale waves has not really been investigated in any of the 
theories, including the wave packet theories of the MJO. 
Here we present a simple MJO model with explicit 
triggered (non QE) convection that emphasizes the 
multi-scale structures of the MJO. 

General circulation models (GCMs) simulate tropical 
intraseasonal variability with varying degrees of fidelity. 



Intercomparison studies (Lin et al. (2006)) show that most 
GCMs underestimate the MJO variance. The weak MJO 
signals in GCMs are believed due to inadequate convection 
schemes in the GCMs. Holloway et al. (2012) compare 



limited-area simulations of the tropical atmosphere over a 
very large domain at different horizontal resolutions with 
both parameterized and explicit convection versions for a 
10 day MJO case study in April 2009. They claim that the 
parameterized models consume CAPE and reach radiative- 
convective equilibrium too quickly, and lack the ability to 
transition from suppressed to active conditions and vice 
versa. As a result, propagating MJO signals are absent in 
the parameterized runs. 

It is widely believed that improvements in the represen- 
tation of subgrid-scale processes in the model would lead 
to a more accurate MJO depiction. However, many of con- 
ventional convection parameterizations are based on the 
quasi-equilibrium (QE) idea, which tends to keep the tem- 
perature profile close to the moist adiabat, and tends to 
damp small-scale waves. A major weakness of parameteri- 
zations is that they artificially separate subgrid-scale pro- 
cesses (departures from QE) from the large scale processes, 
although the small- and large-scale processes interact in na- 
ture. 



waves, which are crucial to the MJO. 

There are two common ways of treating convection; 
triggered and statistical equilibrium convection (or QE). 
There is not a clear distinction between these two cate- 
gories. Over a long period, in the tropics, the generation 
of convective available potential energy (CAPE) by large- 
scale processes nearly balances its consumption by convec- 
tion. The convection may be considered to be in a state 
of statistical equilibrium with the large-scale circulation. 



It is first applied by Arakawa and Schubert (1974). This 



idea has been validated using observational data sets by 



Xu and Emanuel (19891 and Holloway and Neehn (2007). 



QE is a very good assumption to study large-scale circu- 
lations that vary slowly with time compared with convec- 
tive timescales. Such circulations include tropical cyclones. 



the Hadley cell and monsoon circulations {e.g., Emanuel 



et al.| ( |1994[ ), |Emanuel| ( |2007[ )), but so far the MJO has 
not been successfully simulated under the QE context. If 
we are interested in high-frequency small-scale waves, how- 
ever, QE does not work well. First, over a short period, 
CAPE builds up. When some threshold is reached, convec- 
tion is triggered and CAPE is released. Second, in a QE 
scheme, convection will damp small-scale waves faster and 
leave the large-scale waves. We are investigating triggered 
convection in order to simulate high-frequency small-scale 
waves. 

Motivated by the success of SP-CAM simulations, we 
develop a theory emphasizing the role of high-frequency 
small-scale waves. In this paper, we use a shallow wa- 
ter model with triggered convection and radiation repre- 
sented as Newtonian relaxation. Slowly eastward propa- 
gating (MJO-like) signals are observed in our simulations. 
Instead of thinking of the MJO as a large-scale stable or 
unstable mode, we propose the MJO-like signals are inter- 
ference patterns of westward and eastward inertia gravity 
(WIG and EIG) waves. Eastward propagation is due to 
the zonal asymmetry of the EIG and WIG waves. In sec- 
tion 2, we will introduce the shallow water model used in 
this study. In sections 3 and 4 we show our simulation 
results and attempts to understand the simulated signal. 
In section 5, we will discuss our results, and present our 
conclusions and future work. 



Recent studies show improved MJO simulations with 
superparameterizcd (SP) Community Atmosphere Models 
(CAM) (Benedict and Randalll ( |2009l )). The SP-CAM re- 
places conventional boundary layer, moist convection pa- 
rameterizations with a cloud-resolving model embedded in 



each CAM grid ceh ( Khairoutdinov and Randall (2001)) 



The SP-CAM deals with sub-grid scale variability more ac- 
curately. Our interpretation of the SP-CAM results is that 
the cloud-resolving model excites high-frequency small-scale 



2. Model description 

Shallow water model is a two-dimensional model. We 
simulate the upper troposphere by assuming the first baro- 
clinic mode, since the large-scale circulation associated with 
the MJO shows the first baroclinic structure. Thus diver- 
gence in the model refers to upper level divergence and low 
level convergence. Similarly, large layer thickness in the 
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model corresponds to high pressure aloft and low pressure 
near the surface. In this section, we will introduce the shal- 
low water model used in this study. 

This shallow water model describes the evolution of 
constant density, incompressible fluid over the surface of 
the sphere. The model equations are: 



addition to the planetary radius and rotation, there are 
four parameters: the equilibrium geopotianl (/)c, the area 
of one convection event Ag, the convective timescale Tc, 
and the number density of convection events N, which is 
controlled by r through the relation TV ~ . Both Qo and 
r are small, so that only the ratio between the two terms 
matters. 



dtu ^ jv -dxu + -deu H n^xq), 

acost) a a acosa 



(1) 



dtv = -fu ^^dxv + -dev + _ -decf), (2) 

acosa a a a 



dtcp - -V • {V<i>) + s, 



(3) 



(4) 



Equations ([!]) and ([2| are momentum equations, where 
u and V are zonal and meridional velocities; represents 
geopotential corresponding to the layer thickness; a is the 
earth radius; / is the Coriolis parameter (a.k.a. planetary 
vorticity); A and 9 represent longitude and latitude in radi- 
ans. Equation ([3]) is the continuity equation, and s includes 
both mass sources and sinks. In Eq. Q, q represents con- 
vective heating, which is a mass source; and r represents 
radiative cooling, which is a mass sink. 

In this model, convection events are triggered by a low 
value of the layer thickness, i.e. if the layer thickness is 
lower than a threshold 0c, convection will start to add mass 
into this shallow water system: 



10 
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R2) 



when (j) < 
otherwise, 



< 



(5) 



where go is a free parameter of the heating amplitude, and 
Tc is the convective timescale. Each convective event op- 
erates in a certain area = ttR^, where R is the radius 
of each convective event. L = \J Ax^ -I- Ay^ measures the 
distance from the convective center. Radiative cooling is 
constant in both time and space. The ratio |^ determines 
N, the average number of convection events per unit area, 
per unit time. In a statistically steady state, the total mass 
of this system will not change with time. Convection, the 
mass source, will be balanced by the mass sink, radiation. 
The equilibrium geopotential is ^ 0c- 

The SW equations have characteristic length and time 
scales through the planetary radius and the rotation. In 



The Kelvin wave speed, c is equal to v^c- this study, 
we fix c ^ 16 TO s~^ . The mean depth h^. (a.k.a. equivalent 
depth) of this SW system is given by — . Radius of con- 
vection i? is 3 degrees of latitude, whic^ is approximately 
the size of the grid in T42 simulations. Parameter values 
of our control simulation are documented in table 1. We 
vary the horizontal resolution from T42 to T170, and find 
our main results are not sensitive to resolution. The results 
presented in this paper are mainly from T42 simulations. 
We solve these equations in spherical coordinates by using 
the spectral dynamical core of the Geophysical Fluid Dy- 
namics Laboratorys (GFDLs) Flexible Modeling System. 
Although we do not have a moisture variable explicitly in 
our model, our model does illustrate the importance of 
moisture to the MJO. 



3. Simulation results 

Figure [T] shows HovmoUer diagrams of our shallow wa- 
ter simulation from day 500 to day 600. Figure 1 shows 
the symmetric components of the gepotential, zonal wind, 
and convective heating and the antisymmetric components 
of meridional wind. These symmetries and asymmetries 
are respect to the equator. Such meridional symmetry ex- 
cludes the even meridional- wind modes and leaves the odd 
modes. Figure [l^ shows the geopotential. There are two 
major large-scale events labeled A and B. These are the 
MJO-like signals. They move eastward at ~ 3.2 m s~^. 
Sn?all 4©^^ wi^^ are present, and they include Kelvin 
waves, IG waves and Rossby waves. Since the IG waves 
are small scale waves, the absolute values of their speed 
are close to the Kelvin wave speed, which is ^ 16 m s^^. 
Figure[l]D shows the zonal wind. The white represents west- 
ward zonal wind, and the black represents eastward zonal 
wind. The edge between the white and black indicates the 
divergence of zonal wind. Two regions of large scale diver- 
gence are observed, and they are collocated with events A 
and B in Fig. [T^. Figure [TJ: shows the meridional velocity 
V. Small-scale westward and eastward waves are observed, 
but large scale envelopes are not clear in the v field. Figure 
[TJl shows convective heating. Convection is a small-scale, 
short-period process, but two organized long-lasting events 
are observed, and they are collocated with events A and B 
in Fig. [T^ and with large-scale divergence in Fig. [TJd. Sim- 
ilar to the observation, large-scale divergence is collocated 
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with convective centers, and dynamical fields are coupled 
to convection. 

Figure is the HovmoUer diagram of the zonal wind 
from a T42 simulation. Figure 2 shows the HovmoUer di- 
agrams of the zonal wind from T85 and T170 simulations 
with the same parameters. The size of the convection 
in kilometers is the same. They both show propagation 
speeds and horizontal structures that are similar to the 
T42 simulations. This comparison suggests that the T42 
simulations have already converged. In the rest of the pa- 
per, we will present the simulation results from the T42 
simulations. 

To understand the multi-scale structures in our simu- 
lation, we carried out space-time spectral analysis as pio- 



neered by [Wheeler and Kiladis (1999) (hereafter WK99) 



signal. The EOFs together with the corresponding prin- 
cipal components (PCs) show the propagation behavior of 
the MJO-like signal. Combining all the phases of the MJO- 
like signal, we get the horizontal structure shown in Fig. 
[4] The contours represent the geopotential and the vectors 
represent the wind. At the equator, the wind is more zonal, 
and the contours are more parallel to the equator. Away 
from the equator we can see cyclonic and anticyclonic vor- 
tices. In this figure, a wavenumber 2 pattern stands out, 
i.e. the zonal wind alternates from eastward to westward 
twice in the domain. Consider one cycle of the pattern, 
that between 100° and 280° longitude. The maximum di- 
vergence of the zonal wind occurs at 180°. To the west 
of the maximum divergence, at about 160°, there are anti- 
cyclones to the north and the south, collocated with high 
geopotential anomalies. To the east of the divergence of the 
zonal wind, there are cyclones centered at ~ 10° latitude, 
together with low geopotential anomalies. This structure 



is referred to as the quadrupole vortex structure (Majda 



Figures [3^ and b show the power spectra of the symmet- 
ric and anti-symmetric components (about the equator) 
of the zonal wind. Superimposed curves represent disper- 
sion relations of equatorial waves for ~ 16 m Kelvin simulations are confined more closely to the equator than 
wave speed. The dispersion curves of the equatorial waves 
are first derived by Matsuno ( 1966 ) . Different equatorial 



and Stechmann (2009)). The quadrupole vortices in our 



in the MJO as reported by Kiladis et al. ( 2005 1 , where the 



waves are characterized by different dispersion relations 
and meridional mode numbers n. The spectral power under 
the superimposed curves is associated with the correspond- 
ing equatorial waves. In Fig. [3^, we can see spectral peaks 
associated with n = 1 Rossby, Kelvin and n = 1, 3 IG 
waves. Consistent with the slowly eastward moving signals 
in Fig. [ij Fig. |3^ has an MJO-like signal within the white 
box. This signal has planetary scale and low frequency. In 
the rest of this paper, we will try to understand this inter- 
esting phenomena. In Fig. [3jD, we can see spectral peaks 
associated with n = 2 Rossby, MRG and n = 2 IG waves. 
A striking feature of Fig. [3] is the intense power associ- 
ated with high frequency IG waves. High IG wave activity 
is associated with the convective parameterization in our 
model. We will argue that high IG wave activity is the 
key to the MJO-like signal. Most of the power is concen- 
trated in low wavenumber and low frequency region, except 
for the spectral peaks associated with high frequency IG 
waves. This is, in general, a red noise spectrum. Although 
different in details. Fig. [3] captures some fundamental fea- 
tures of the observed spectra by WK99. 

Figure |4] shows the horizontal structure of the MJO- 
like signal in our simulation. To get Fig. |4j we first get 
the MJO filtered signal through zonal wind field in the 
wavenumber-frequency domain. The filtering excludes all 
wavenumbers and frequencies except those in the MJO box 
(Fig. |3^). Following Wheeler and Hendon (2004), we car- 



ried out empirical orthogonal function (EOF) analysis of 
the MJO filtered signal. We found the first two EOFs can 
contribute ~ 85% of the total variance of the MJO filtered 



off-equatorial vortices center at ~ 20° latitude. 



4. Proposed mechanism 

MJO-like signals have been simulated in our model. 
The simulated signal captures major features of the MJO, 
including propagation speed, horizontal scale, multi-scale 
structures and quadrupole vortex structures. Then we 
need to understand what the energy source is to main- 
tain the long-lasting signal, and also need to address the 
eastward propagation mechanism. Figure [5] shows high res- 
olution HovmoUer diagrams. In this high resolution view, 
we are able to diagnose how convection is triggered, and 
how waves are excited. In Fig. [5^, there is a standing oscil- 
lation (SO) at a longitude of - 206°. On day 163.75, there 
is a local minimum of at 206°. This triggers a convec- 
tion event. A quarter cycle later, the convection reaches its 
maximum value, and a quarter cycle after that, phi reaches 
a local maximum. Such cycles repeat 5 times in Fig. [5]with 
a period of 0.5 day, which is twice r^. The amplitude of 
convective heating is small, so there is no nonlinearity due 
to advection. Linear waves are excited. However, due to 
the positive-only heating, nonlinearity has been introduced 
into this system. Net zonal divergence of u collocates with 
the SO event. The SO is a representation of selectively am- 
plified waves. Thus energy is introduced into this shallow 
water system. The SO can also be viewed as a stand- 
ing wave, composed of interfering eastward and westward 
waves with similar speeds. In the power spectra, WIG and 
EIG waves are the only waves that propagate toward each 
other with same meridional structures and similar prop- 
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agation speed. Therefore, we propose that the simulated 
MJO signal is an interfering pattern between WIG and 
EIG waves, which are excited by convection. The propaga- 
tion speed of the MJO pattern is associated with the phase 
speed difference between WIG and EIG waves. 

Figure [6] shows how the simulation results matches the 
theoretically derived propagation speed. If the hypothesis 
is correct, the MJO propagation speed is one half of the 
phase speed difference between the WIG and EIG waves. 
In Eq. ([6| , the EIG phase speed is c\ , the WIG phase speed 
is C2, and the cosine factor on the right is the standing 
oscillation. The sine factor on the right is the drift of the 
pattern at the MJO propagation speed. 

sm[fc(a; — cit)] + sin\k(x + C2t)\ = 2sin[k(x t)]cos 

^ (6) 

The thick black lines denote the theoretically derived MJO 
propagation speed from the dispersion relation of IG waves, 
and the markers represent the results from simulations with 
different convective timescale. The lower thick line is for 
the first meridional mode, and the higher thick line is for 
the third meridional mode. The abscissa is the correspond- 
ing frequency of IG waves. The theoretically derived speed 
does not fit the simulated results perfectly, especially for 
short convective timescales (converted into frequency) . In- 
creasing the convective timescale, the MJO propagation 
speed increases monotonically, and this behavior is qualita- 
tively consistent with the theoretically predicted behavior. 
The marker starts to locate between the two thick lines 
when Tc is longer than 0.6 day, i.e., the simulated MJO 
signal propagates faster than the predicted speed from the 
first meridional mode, but slower than that of the third 
meridional mode. This suggests either that convection am- 
plifies higher meridional modes at higher frequencies or 
that convection excites lower frequencies than It is 

possible that due to the positive-only convective heating, 
waves with frequency lower than tt— will be amplified. 

As an interference pattern, the horizontal scale of the 
MJO will be associated with the bandwidth of the excited 
waves. Our simulations show that the bandwidth of waves 
is associated with how frequently the convection happens. 
If convection is very frequent, a wide band of waves will be 
excited. However, if there are only a few convection events 
during the simulation, only a narrow band of waves will be 
excited. Since the convection frequency per unit area N is 
of order decreasing r will decrease N if all the other 
parameters are fixed. As a result, the MJO wavenumber 
decreases. Figure 7 shows the results of systematic ex- 
periments where the MJO scale increases as r increases. 
The amplitude of the convective forcing is small so that 



the waves are linear waves, and increasing the convective 
strength does the same job as decreasing r. 

5. Discussion and Conclusions 

In this paper, we have presented a shallow water model 
with triggered convection and simple radiation treatments. 
MJO-like signals are observed in our simulations. Our sim- 
ulation results are robust over a wide parameter range, and 
do not depend on resolutions. We further propose that the 
MJO-like signal is an interference pattern of the WIG and 
EIG waves, whose frequency is set by the duration of in- 
dividual convective events. We have tested this hypothesis 
by systematically varying the parameters Tc and N. Both 
the propagation speed and the zonal wavenumber of the 
J]Oey}fc signals respond to the parameter changes con- 
sistgntly -with the theoretical predictions. Instead of the 
structure of the IG waves, the MJO-like signal exhibits the 
quadrupole vortex structures. This might indicate that 
these off-equatorial vortices are responses to the large-scale 
envelope of convective heating near the equator, rather 
than free wave structures. Our simulation results suggest 
that the MJO perhaps is not a large-scale low-frequency 
wave in which, convection acts as QE adjustment. Small- 
scale high-frequency waves might be crucial. Therefore, in 
order to simulate the MJO, the behavior of moist convec- 
tion in a short time period should be represented properly. 

We use a one-layer shallow water model without explicit 
moisture variables. However, it does not mean that we ex- 
clude the effect of moisture for the following reasons. First, 
triggered convection only occurs in the moist atmosphere, 
where conditional instability could exist. In a shallow wa- 
ter model, the geopotential (j) is a measure of the static 
stability N'^ of the atmosphere. Low indicates the re- 
duced iV^. Only if iV^ has been reduced to a critical value, 
in other words, CAPE has been accumulated by a certain 
amount, convection will be triggered. This mimics the pro- 
cess of conditional instability and triggered convection in a 
moist atmosphere with the minimum recipes. Second, we 
set the Kelvin wave speed close to the convectively coupled 
Kelvin wave speed, which is '-^ 16 to s~^. 

Due to the simplicity of our model, we have to make as- 
sumptions about the MJO vertical structure and the Kelvin 
wave speed. Although these assumptions are consistent 
with observations, a complete theory will have to explain 
why the MJO has a first baroclinic structure, and why the 
Kelvin wave speed is ~ 16 to s^^. The next step is to 
use a 3D model to test our hypothesis. We do not need 
to assume the MJO vertical structure and the Kelvin wave 
speed in the 3D model. To keep the key features of our 
current model, we will implement triggered convection in 
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the 3D model. In this model, CAPE will be accumulated 
until convective inhibition is close to zero. Then CAPE 
will be released. 

It is also worth to look for IG waves and their rela- 
tions to convection and the MJO in observations to test 
our hypothesis. The correlation between the MJO and 
the IG waves has already been examined using observa- 



tional datasets {e.g. Yang and IngersoU (20111, Yasunaga 



and Mapes (2011)). Statistically significant correlations 
have been identified, but the IG waves can only explain 
very limited variance of the MJO. However, previous stud- 
ies do not help to test our hypothesis. First, most of the 
previous studies used coarse temporal and spacial resolu- 
tion datasets, so they cannot resolve high-frequency IG 
waves. Second, previous studies examined the correlations 
between the MJO and the WIG waves, and the MJO and 
the EIG waves separately. The interference pattern of the 
WIG and EIG waves are not even included. To test our hy- 
pothesis, one should first combine the WIG and EIG signals 
in a high-resolution dataset, and identify their interference 
patterns. Then examine the correlations between the MJO 
and this identified interference pattern. 
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Table 1. Parameter values in the control simulation. 
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Fig. 1. HovmoUer diagrams of the (a) geopotential, (b) zonal wind, (c) meridional wind and (d) convective heating of our 
shallow water simulation from -15° to 15° latitude. The geopotential, zonal wind and convective heating are symmetric 
components about the equator, and the meridional wind is the anti-symmetric component. The white represents low, 
and the black represents high. 
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Fig. 2. Hovmoller diagrams of zonal wind from (a) T85 simulation, (b) T170 simulation. The white represents low, and 
the black represents high. 



11 




Zonal Wavenumber 




-40 -30 -20-10 10 20 30 40 



Zonal Wavenumber 



Fig. 3. Zonal wavenumberfrequency power spectra of zonal wind from -15° to 15° latitude, (a) is for the symmetric 
component, and (b) is for the anti-symmetric component. Red represents high power density, and blue represents low 
power density. Red, blue, black and green lines denote dispersion curves of IG, Rossby, Kelvin and MRG waves for 
different meridional modes. The white box ranges from 1 to 5 in wavenumber, and from ^ to cycle per day (CPD) in 
frequency. 
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Fig. 4. Longitude-latitude map of eomposites of tfie MJO-like signals. The arrows denote the wind field V = {u, v) , and 
the contours denote the geopotential. The thick solid (thin dashed) contours represent positive (negative) anomalies. 
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Fig. 5. HovmoUer diagrams of the (a) geopotential, (b) zonal wind, (c) meridional wind and (d) convective heating of our 
shallow water simulation from -15° to 15° latitude. The geopotential, zonal wind and convective heating are symmetric 
components about the equator, and the meridional wind is the anti-symmetric component. The blue contours denote 
low, and the red contours denote high. 
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Fig. 6. Propagation speed (m s~^) versus frequency (CPD). The dark lines are derived from the dispersion relation of 
WIG and EIG waves. The lower one is for the first meridional mode, and the upper one is for the third meridional mode. 
The markers represent the simulation results for different convective timescales (tc). We convert to frequency by using 
frequency = . 
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Fig. 7. Zonal wavenumber versus normalized number density of convection where N — No is given in Table jlj 
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